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We study the fluid squeeze-out from the interface between an elastic solid with a flat surface and 
a rigid solid with a randomly rough surface. As an application we discuss fluid squeeze-out between 
a tire tread block and a road surface. Some implications for the leakage of seals are discussed, and 
experimental data are presented to test the theory. 



1. Introduction 

Contact mechanics between solid surfaces is the ba- 
sis for understanding many tribology processes [ll-H] such 
as friction, adhesion, wear and sealing. The two most 
important properties in contact mechanics are the area 
of real contact and the interfacial separation between 
the soUd surfaces. For non-adhesive contact and small 
squeezing pressure, the average interfacial separation de- 
pends logarithmically on the squeezing pressure [1, 
and the (projected) contact area depends linearly on the 
squeezing pressure [lO|. Here we study how the (average) 
interfacial separation depends on time when two elastic 
solids with rough surfaces are squeezed together in a fluid. 
In particular, we calculate the time necessary to squeeze- 
out the fluid from the contact regions between the solids. 
As an application we discuss fluid squeeze-out between a 
tire tread block and a road surface. Some implications 
for the leakage of seals are discussed, and experimental 
data are presented to test the theory. 

2. Squeeze-out: large separation 
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FIG. 1: The upper solid, a cylindrical block with the diam- 
eter D = 2R, the elastic modulus E and the Poisson ratio 
V, is squeezed in a fluid against a rigid substrate. The bot- 
tom surface of the cylinder is perfectly flat and the substrate 
surface randomly rough with the root-mean-square roughness 
amplitude /irma- 



Consider an elastic solid with a flat surface squeezed 
in a fluid against a rigid solid with a randomly rough 
surface, see Fig. [T] The fluid is assumed to be Newtonian 
with the viscosity r/. The upper solid is a cylindrical 
block with the radius R, the elastic modulus E and the 
Poisson ratio v. The bottom surface of the cylinder is 
perfectly flat, and the substrate randomly rough with 
the root-mean-square roughness amplitude h^ms- Here 
we focus first on the the simplest possible situation which 
can be studied analytically, where the (macroscopic or 
locally averaged) pressure distribution in the fluid gives 
rise to negligible deformations of the bottom surface of 
the elastic block. This requires that the amplitude 5u of 
the (fluid-induced) elastic deformations is much smaller 
than /iims- Since Su is typically of order (or smaller than) 
~PqR/E*, where E* — E/{\ — i''^) andpo is the pressure 
applied to the upper surface of the cylinder block, we 
get the condition h^^s » PoR/E* . For elastically stiff 
materials with E* « 10^^ Pa and for R = \ cm we get 
PoR/E* < 0.1 /im if Po < 1 MPa. For tread rubber E* « 
10 MPa and if po « 0.3 MPa we get poR/E* « 0.3 mm 
which is smaller than the root-mean-square roughness 
of many asphalt road surfaces. In many applications a 
thin rubber film, coating a hard solid, is in contact with 
an elastically hard countersurface. If the linear size of 
the contact region is large compared to the rubber film 
thickness, this geometry will strongly suppress the (fluid- 
induced) deformations of the rubber film on the length 
scale of the linear size of the nominal contact area, and 
most of the (non-uniform) deformations of the rubber 
film is due to the interaction with the substrate asperities. 

We first develop a theory which should be accurate for 
large enough interfacial separation, e.g., corresponding to 
the early phase of the squeeze-out process. We assume 
that the longest wavelength roughness component, Aq, 
is small compared to the linear size R of the (apparent) 
contact region. In this case we can speak about locally 
averaged (over surface areas with linear dimension of or- 
der Ao) quantities. 

Neglecting inertia effects, the squeeze-out is deter- 
mined by (see, e.g., Ref. Q) 

du ^ 2fl^pfluid(^) .-.s 
dt ~ 3r/i?2 ' ^ > 

where p^u\d{t) is the (average) fluid pressure, and u the 



(locally averaged) interfacial separation. If po is the ap- 
plied pressure acting on the top surface of the cylinder 
block, we have 



Pfluid(0 ^ Po- Pcont(t), 



(2) 



where Pcont is the asperity contact pressure. We will 
first assume that the pressure po is so small that for all 
times u » /irms- In this case we can use the asymptotic 
relation 81 



(3) 



Pcont ~ PE exp . 

"0/ 



where uq = /irms/ck. The parameters a and^ depends 
on the fractal properties of the rough surface[8|. 
From (3) we get 



du 
It 



Uq dpc 



Pcont dt 

Using (4) and (2) we get from (1): 



dt 



Pcont {PO - Pcont) 



(4) 



(5) 



For long times Pcont ~ Po and we can approximate (5) 
with 



dpcont (po) 



dt 3r/R'^uo 
Integrating this equation gives 

Pcont (i) ~ Po - bo - Pcont (0)] CXp 



Po(PO - Pcont) 




where 



2/i?ms-Po 2aft2^,po ■ 



(6) 



Using (3) this gives 



Pcont (0) 



Po 



uoexp 




where Uoc = Uf)\og{l3 E* / po) . Thus, u{t) will approach 
the equilibrium separation Uoo in an exponential way, 
and we can define the squeeze-out time as the time to 
reach, say, l.Oluoo- For flat surfaces, within continuum 
mechanics, the film thickness approach zero as t — ?> cxd as 
u ^ i"^/^. Thus in this case it is not possible to define a 
meaningful fluid squeeze-out time. 

Let us measure distance u in units of /irms, pressure in 
units of Po and time in units of r (eq. (6)). In these units 
(5) takes the form 



dpcont _ , -3 



dt 



U Pcont (1 - Pcont) , 



(7) 




FIG. 2: The interfacial separation (or film thickness) u (in 
units of ftrms) as a function of the squeeze-time t (in units of 
r) for several values of the parameter k. For each squeeze- 
pressure, the upper (red) lines are the result for the rough 
surfaces while the lower (blue) lines are for flat surfaces. 



In the same units (3) takes the form 

Pcont ~ K~^exp {-au) . 



where 



Po 

PE* 



(8) 



(9) 



In Fig. [2] we show the interfacial separation (or film 
thickness) u (in units of /irms) as a function of the 
squeeze-time t (in units of t) for several values of the 
parameter n. For each k value, the upper (red) lines are 
the result for the rough surfaces while the lower (blue) 
lines are for fiat surfaces. In the calculation we have used 
a = 0.5 and P — 1. 

In Fig. |3] we show the fluid squeeze-out time t* (in 
units of r) and the final interfacial separation (or film 
thickness), u{t*) (in units of /irms), as a function of the 
parameter k. We define t* so that u{t*) = l.Ofu(oo). 

At high enough squeezing pressures, the interfacial sep- 
aration after long enough times will be smaller than /irms, 
and the asymptotic relation (3) will no longer hold. In 
this case the relation Pcontiu) can be calculated using the 
equations given in Ref. il[ (see also Appendix A). Sub- 
stituting (2) in (I) and measuring pressure in units of 
Po, separation in units of /ir 
obtain 



and time in units of r one 



du 
'di 



'^U^(I - Pcont), 



(10) 



where a = /irms/uo- This equation together with the re- 
lation Pcont (w) constitute two equations for two unknown 
(m and Pcont) which are easily solved by numerical inte- 
gration. In what follows we refer to the theory presented 
above as the average-separation theory. 

3. Squeeze-out: general theory 




FIG. 3: The fluid squeeze-out time t* (in units of r) and the 
final interfacial separation (or film thickness), u{t*) (in units 
of /irms), as a function of the logarithm (with 10 as basis) of 
the parameter k. We define f so that u{t*) — l.Olit(oo). 



We now present a general theory of squeeze-out, which 
is accurate for small separation and which reduces to the 
result presented in Sec. 2 for large separation. The the- 
ory presented below is based on a recently developed the- 

We assume again 



ory of the leak-rate of (static) seals [11 



that the longest wavelength roughness component, Aq, is 
small compared to the linear size R of the (apparent) 
contact region. In this case we can speak about locally 
averaged (over surface areas with linear dimension of or- 
der Ao) quantities. Let J(x, t) be the (locally averaged) 
2D-fluid flow vector which satisfies the continuity equa- 
tion 



V- J 



du 
'dt 



0, 



(11) 



where u(x, t) is the (locally averaged) surface separation 
or, equivalently, the 2D-fluid density (fiuid volume per 
unit area). Here and in what follows 



{Jx: Jy 



are 2D vectors. In Ref. 



(12) 



11 



we have shown that within 



an effective medium approach 



-CTcffVpflu 



(13) 



where Pfluid(x, is the (locally averaged) fiuid pressure 
and where the effective conductivity croff(Pcont) depends 
on the (locally averaged) contact pressure Pcont(x, i). 
Note that when inertia effects are negligible 



d?x [pcont(x,t) -KMuid(x,i)] = F^{t) (14) 



is the applied normal load. The function crcff(pcont) can 
be calculated from the surface roughness power spectrum 
and the (effective) elastic modulus as described in Ref. 



ll|. Substituting (12) and (13) in (11) gives 

du 
'dt' 



uidj 



(15) 



Eq. (15) together with the relation Pcont('S) and the 
(standard) expression relating the macroscopic deforma- 
tion u{x.,t) to the local pressure po(x) = Pcont(x,t) + 
Pfluid(x, i) constitutes three equations for the three un- 
known Pcont, Pfluid and u. In addition one need the effec- 
tive medium expression for tToft (pcont ) , and the "bound- 
ary condition" (14) must be satisfied. Here we will not 
study the most general problem but we focus on the lim- 
iting case discussed above where the macroscopic defor- 
mations of the solid walls can be neglected. In this case u 
and pcont will only depend on time. As a result a-eff(Pcont) 
will only depend on time. Thus, (15) reduces to 



du 
'dt' 



Since the right hand side only depends on time. 



Pfluid(x,t) = 2pfluid(i) 1 - {r/R) 



(16) 



(17) 



where Pfiuid(i) is the average (nominal) fiuid pressure in 
the nominal contact region. Substituting (17) in (16) 
gives 



'CTeff 



i?2 

Eq. (14) takes the form 



-Pfluid 



du 
~di' 



Pcont + Pfluid = PO 



Using (4) and (19) in (18) gives 



dpc 



8CToff(Pc 



dt 



R^uo 



-Pcont(PO -Pcont). 



(18) 



(19) 



(20) 



From this equation one obtain Pcont{t), and using (3) 
and (19) one can then calculate u(t) and pfluid(i)- As 
shown in Appendix A, when pcont — > 0, w — )■ oo and 
(Tcff — > u? /\2rj, which is an exact result to leading order 
in hj-ms/u. Substituting acS — u^/12ri in (20) gives (5). 
Thus, in the limit of small pressures po, the present treat- 
ment reduces to the average-separation theory of Sec. 2, 
which is exact when the average separation u between 
the surfaces is large (which is the case for all times if the 
pressures po is small). We will refer to cTcs = /12rj as 
the average-separation expression for acS- 

Let us study the squeeze-out for long times. For long 
times Pcont ~ Po and we can approximate (20) with 



dpc 



^o-off(Po 



Po(PO -Pcont). 



dt R'^uq 
Integrating this equation gives 

Pcont{t) ~ Po - [PO - Pcont(O)] CXp ( - 



(21) 



8 gcfF(Po)PO _ 

R^uo 



so that Pcont(i) approach po (and u{t) approach Uoo) in 
an exponential way, just as for the simpler model studied 
in Sec. 2 

If we measure pressure in units of po, separation u in 
units of hrms, and time in units of t, Eq. (20) takes the 
form 



dpc 



dt 



O" offPcont ( 1 — Pcont ) , 



where 



aoff = 12-qacs/h 



3 

rins ■ 



(22) 



(23) 



The relation between pcont and u is given by (3). 

At high enough squeezing pressures, the interfacial sep- 
aration after long enough times will be smaller than /irms, 
and the asymptotic relation (3) no longer hold. In this 
case the relation Pcont(w) can be calculated using the 
equations given in Ref. [1] (see also Appendix A). Sub- 
stituting (19) in (18) and measuring pressure in units of 
Poj separation in units of /irms and time in units of r one 
obtain 



du 



Vcff (1 — Pcont) 



(24) 



This equation together with the relations Pcont('«) and 
Ceff (Pcont ) constitute three equations for three unknown 
{u, Pcont and (Teff) which are easily solved by numerical 
integration. In the critical junction theory which we will 
use below CTcff = (a"i(Cc))f /12?7, where the separation 
Mi(Cc) in defined in Ref. [11| (see also Appendix A) and 
where a < 1 is a number of order unity (see Ref. [11]). 

Fig. m shows the calculated interfacial separation u{t) 
as a function of the squeeze-time t for a silicon rubber 
block (elastic modulus E = 2.3 MPa) squeezed against 
a rough copper surface (log- log scale with 10 as basis), 
with the power spectrum given in Ref. [l2l |. Curves 1 
and 2 are the theory predictions using (10) and (24), re- 
spectively. In (24) we have used cTeff as calculated using 
the critical-junction theory described in Ref. (ll| . which 
gives nearly the same result as the effective medium the- 
ory described in the same reference. Results are also 
shown for a flat substrate surface. The rubber block is 
assumed to be cylindrical with the radius i? = 1.5 cm 
and the surface of the copper block has the root-mean- 
square roughness 0.12 mm. The squeezing pressure 
Po = 0.2 MPa (upper curves) and 1.4 MPa (lower curves), 
corresponding io k — po/f3E* w 1.5 and 10.3, respec- 
tively. Note that for the pressure po — 1.4 MPa, after 
long enough time the area of real contact. A, percolate 
(i.e., A/Aq > 0.5, see Ref. Um), and there is no fluid 
leak channel at the interface [13j. As a result, when the 
contact area percolate the fluid is confined at the inter- 
face and is not able to leak-out. Thus, even after very 
long time the interfacial separation is larger than would 
be expected in the absence of trapped or confined fluid 



J pressure 
/ 0.2 MPa 




FIG. 4: The interfacial separation u{t) as a function of the 
squeeze-time t for rough copper surface (log-log scale with 
10 as basis). Curves 1 and 2 are the predictions using the 
average-separation theory and using the critical-junction the- 
ory, respectively. Results are also shown for a flat surface. 
The rubber block is assumed to be cylindrical with the ra- 
dius R — 1.5 cm. For a copper surface with the root- 
mean-square roughness 0.12 mm. The squeezing pressure 
Po = 0.2 MPa (upper curves) and 1.4 MPa. The elastic mod- 
ulus E = 2.3 MPa. 
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FIG. 5: The normalized contact prsssur© pcon (t) /po as a func- 
tion of the logarithm (with 10 as basis) of the squeeze-time t 
for rough copper surface. Curves 1 and 2 are the predictions 
using the average-separation and using the critical-junction 
theory, respectively. The rubber block is assumed to be cylin- 
drical with the radius R = 1.5 cm. For a copper surface 
with the root-mean-square roughness 0.12 mm. The squeez- 
ing pressure po = 0.2 MPa (upper curves) and 1.4 MPa. The 
elastic modulus E = 2.3 MPa. 



(e.g., for dry contact), where no part of the load would 
be carried by the fluid. 

Fig. [5] shows the normalized contact pressure 
Pcon{t)/po as a function of the logarithm of the squeeze- 
time t for the same system as in Fig. Curves 1 and 
2 are the predictions using (10) and (24), respectively. 
Note that for the pressure po = 1.4 MPa, even after very 



long enough time Pcon < 0.9po- This is again a conse- 
quence of the fact that the non-contact area does not 
percolate and fluid is confined at the interface, and even 
after very long time, more than 10% of the external load 
is carried by the confined fluid. 

Following Ref. the analysis presented above may 
be extended to include the fluid-pressure induced elas- 
tic deformation of the solid surfaces at the interface. It 
is also easy to include the dependency of the fluid vis- 
cosity on the local pressure or local surface separation. 
The former is important for elastically hard solids (high 
pressures), e.g., steel, and the latter "confinement effect" 
even for elastically soft solids (low pressures) when the 
fluid film thickness becomes of order a few nanometer or 

Finally, let us note the following: The interfacial sep- 
aration is usually mainly determined by the most long- 
wavelength surface roughness, which is observed close to 
the lowest magnification C ~ 1; the shortest wavelength 
(small amplitude) roughness has almost no influence on 
u. However, it is also of great interest, e.g., in the con- 
text of tire friction on wet road surface (see Sec. 4), to 
study how the fluid is squeezed out from the (apparent) 
asperity contact regions observed at higher magnifica- 
tion (, see Fig. [T] The theory developed above can be 
applied to this case too. Thus, let us study the squeeze- 
out of fiuid from the apparent asperity contact regions 
observed at the magnification (. At this magnification 
no surface roughness with wavelength below L/^ can be 
observed. However, when the magnification is increased 
one observe shorter wavelength roughness which will in- 
fluence the fluid squeeze-out, and which may even result 
in sealed-off, trapped fluid. We can apply the theory 
above to study the squeeze out of fluid from the asperity 
contact regions observed at the magnification by using 
instead of the external pressure poi the local squeezing 
pressure p{() = poAo/A{(), where A{() is the (appar- 
ent) contact area observed at the magnification ^. The 
surface roughness in the contact regions is given by the 
surface power spectrum C{q) for q > (/L. With these 
modifications we can use the theory above to calculate 
the squeeze-out of fluid from the apparent asperity con- 
tact regions observed at the magnification (. One com- 
plication is, however, that the fluid is squeezed-out from 
an asperity contact region into the surrounding, and the 
fluid pressure in the surrounding may be higher than the 
external pressure (which we have taken as our reference 
pressure in the study above) existing outside the nomi- 
nal contact region observed at the lowest magnification 
C = 1. As a result, in order to study the squeeze-out 
from the asperity contact regions observed at the magni- 
fication C, one must first study the squeeze-out from the 
asperity contact regions observed at lower magnification 
1 <(<('. We will not develop this theory here, but we 
believe a similar approach as that used to describe mixed 
lubrication for flat on flat in Ref. may be applied to 
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FIG. 6: The pressure flow factor (pp, calculated using the 
Bruggeman efi'ective medium theory (see Ref. [lil), as a func- 
tion of the average surface separation u for a rough copper 
surface with the root-mean-square roughness 0.12 mm. The 
squeezing pressure po — 1.0 MPa and the elastic modulus 
E = 2.3 MPa. 



the present problem. 

4. Pressure flow factor 

The theory presented above can be used to calculate 
the pressure fiow factor (pp first introduced by Patir and 
Cheng This quantity is defined so that the (locally 
averaged) flow current J associated with fluid flow at 
the interface between two stationary solids with rough 
surfaces, is given by 



1 

i 

I2ri 



uid 



where u and pauid are the locally averaged interfacial sep- 
aration and the locally averaged (nominal) fluid pressure, 
respectively. Using (13) and (23) we get 

(/>p = 1277Crcff/u^ = CToff(/lrms/'S)^ 

In Fig. iniwe show the pressure flow factor 0p, calculated 
using the Bruggeman effective medium theory (see Ref. 
11| ) , as a function of the average surface separation u for 
a rough copper surface with the root-mean-square rough- 
ness 0.12 mm. The squeezing pressure p^) — 1 MPa and 
the elastic modulus E — 2.3 MPa. Note that (pp vanish 
already for a non-zero u. This is due to the percolation 
of the contact area. Note also that (/)p — >■ 1 as m — oo. 
This result is expected because for large separation the 
surface roughness should have a negligible influence on 
the fluid flow. 

The pressure flow factor (f>p is usually determined by 
calculating (numerically) the fluid flow in small interfa- 
cial units jl6j|. Most studies neglect the elastic deforma- 
tions of the solid walls, and only take into account sur- 
face roughness over two decades (or less) in length scale. 
The present treatment includes elastic deformation and 
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FIG. 7: A rubber block squeezed against a randomly rough 
substrate in a fluid. The (average) interfacial separation u(t) 
decreases with time t due to squeeze-out of the fluid. The 
time dependence of it depends mainly on the long- wavelength 
roughness components. Nevertheless, the short wavelength 
roughness components will affect the squeeze-out in asperity 
contact regions (observed at the magnification which is rel- 
evant for, e.g., rubber friction. This can be studied using the 
same theory applied to an asperity contact region by replac- 
ing the external pressure po with the local squeezing pressure 
p{() at the asperity (see text for details). 



can easily take into account roughness on arbitrary many 
decades in length scale. In Ref. [l3| we have shown how 
one can generalize the treatment above to obtain the flow 
factor (now a matrix) for surfaces with random roughness 
with anisotropic statistical properties. 

5. Application to tire on wet road 

As an apphcation, consider a tire tread block squeezed 
against asphalt road surfaces in water. In Fig. |S] and 
[9] we show the time dependence of the interfacial sepa- 
ration u, and the difference Au = u(t) — u{oo), respec- 
tively. In the calculation we have used the theory of Sec. 
2, which is valid in the present case where (for all times) 
u > 2hrms (see below). We show results for two road 
surfaces, with the root-mean-square roughness 0.72 mm 
(surface 1) and 0.24 mm (surface 2). We assumed the 
squeezing pressure po = 0.2 MPa and the elastic modulus 
_E = 10 MPa. In the calculation we have used the surface 
roughness power spectrum obtained from the measured 
surface topographies. Note that for the smoother sur- 
face the squeeze-out time is roughly one decade longer 
than for the rougher asphalt road surface. In order for 
the water to have a negligible influence on the hysteresis 
contribution to the friction, the water layer in the road- 
rubber contact regions must be smaller than ~ 1 /im. For 



the smoother road surface 2 it takes about '-^3xl0~^s 
to reach Au = 1 /xm. If the tire rolling velocity is 30 m/s 
and the length of the tire foot-print 0.1 m, then a tread 
block spend about 3 x 10~^ s in contact with the road. 
Thus, from the calculation above one may conclude that 
accounting just for the viscosity of the fluid (water) one 
expect (during rolling) almost complete fluid squeeze-out 
from the tread-block road contact area, during most of 
the time the tread block spends in the footprint. How- 
ever, at the start of raining after a long time of dry road 
condition, the water will be mixed with contamination 
particles (e.g., small rubber and road wear particles), and 
the effective viscosity of the mixture may be much larger 
than for pure water. In this case the squeeze-out may be 
incomplete, which could result in viscous hydroplaning 
during braking. In addition, even for pure water, regions 
of sealed off (trapped) fluid may appear at the interface 
at high enough magnification, which will reduce the hys- 
teresis contribution to the tire- road friction[l^. We note 
that during braking at small slip (below the maximum in 
the ^-slip curve) the tread block does not slip until close 
to the exit of the tire-road footprint, and the discussion 
above should therefore be valid for this case too. 

It is interesting to compare the result above with the 
squeeze-out time due to inertia (but neglecting the vis- 
cosity) . The time-dependence of the squeeze-out for flat 
surfaces is given by (see Ref. fl9| ) 



where 



u{t) = u(0)exp(-</T'), 



64pc 



Thus, the time it takes to reduce the film thickness from 
m(0) to a thickness of order h^rna is 

r «r'log[u(0)//i„„s]. 

If u(0) = 1 cm and /irms — 0.3 mm we get t* w 
T'log(30) « 3.4r'. With D « 2 cm, po « 0.3 MPa 
and p w lO'^ kg/m^ we get r' « 1.4 x 10~* s and 
r w 5 X 10--* s. 

6. Leakage of seals 

In the calculation of the leak-rate of seals presented in 
Ref. JTj we neglected the influence of the fluid pressure 
on the contact mechanics. This is a good approximation 
as long as the squeezing pressure po is much higher than 
the fluid pressure pfluid, which was the case in the experi- 
ments presented in Ref. ll|. However, in many practical 



situations it is not a good approximation to neglect the 
influence of the fluid pressure on the contact mechanics. 
Since the fluid pressure is higher on the fluid entrance 
side than on the fluid exit side, one expect the elastic 
wall to deform and tilt relative to the average substrate 
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FIG. 8: The interfacial separation u(t) as a function of the 
squeeze-time t for two asphah road surfaces (log-log scale with 
10 as basis) . The tread block is assumed to be cylindrical with 
the radius R — 2 cm. For two asphalt road surfaces with the 
root-mean-square roughness 0.72 mm (surface 1) and 0.24 mm 
(surface 2). The nominal squeezing pressure po = 0.2 MPa 
and the elastic modulus E = 10 MPa. The vertical lines 
denote (from left to right) the time when u{t)/u{oo) — 1.1, 
1.01 and 1.001. The thin lines denote the fluid film thickness 
for perfectly fiat surfaces. The calculations use the average- 
separation expression for acs- 




FIG. 9: The difference Ait = u{t) — u(oo) as a function of the 
squeeze-time t for two asphalt road surfaces (log-log scale with 
10 as basis) . The tread block is assumed to be cylindrical with 
the radius R — 2 cm. For two asphalt road surfaces with the 
root-mean-square roughness 0.72 mm (surface 1) and 0.24 mm 
(surface 2). The squeezing pressure po = 0.2 MPa and the 
elastic modulus i? = 10 MPa. 



surface plane, see Fig. [TUl Here we show how one can in- 
clude the fluid pressure when calculating the leak rate of 
seals. For simplicity we focus on the simplest case where 
the fluid pressure only depends on one coordinate X, as 
would be the case for most seal applications, e.g., rubber 
0-ring seals, see Fig. [TUl In this case, for a stationary 
situation (15) takes the form 

^ fo-ofr(Pcont(a;))^Pfluid(a;) j = 



o-offbcontla;))— pfluid(a;) = B 
where 5 is a constant. From this equation we get 

Pduidix) =A + b[ dx' a-f^{pcontix')) (25) 
Jo 

where A is a constant. If the fluid pressure for a; < 
(high pressure side) is denoted by pa, and for a; > L (low 
pressure side) by pb, then using that Pfluid(O) — Pa. and 
Pfluid(-^) — Ph we can determine the constants A and B 
in (25) and get A = and 



B = 



Ph -Pa 



lo c^off (Pco„t(a:')) 
Substituting these results in (25) gives 

dx' a-g{p^orA{x')) 



Pfluid(a;) = pa + (pb - Pa) 



/o f^^' f^off (Pcont(a;')) 



(26) 



(27) 



The rubber 0-ring is squeezed against the substrate 
by the normal force per unit radial length, /n, see Fig. 
[TUl In the contact region between the cylinder and the 
substrate occur a nominal (locally averaged) pressure: 



Po{x) = _Pcont(2^) +Pfluid(2:)- 

We consider a stationary case so that 



dx Po{x) 



h 
L ' 



The elastic deformation field 

.2 2 



X 



U{x) = Uc+^ 



nE* 



dx' po{x')log 



(28) 



(29) 



(30) 



Equations (27), (28) and (30), together with the equation 
determining the relation between pcont and u, represent 
4 equations for the 4 unknown variables po, Pcont, Pfluid 
and u. In addition the pressure po{x) must satisfy the 
normalization condition (29) which determines the pa- 
rameter Uc in (30). 

The leak-rate of the seal is given by Q = JxLy where 
Ly is width of the seal (e.g., the circumstance of the seal 
for a rubber 0-ring). Using (13) we get 

Q = -LyC7csiPcont{x))-^Pflnidix) = -LyB 



Using (26) this gives 



LyiPi,-Ph) 



If Pcont is constant this gives 



Q = 7^croff(Pa -Pb) 



(31) 



(32) 
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FIG. 10: Cross-section of rubber 0-ring squeezed against a 
rigid, randomly rough surface in a fluid. The fluid pressure is 
higher at a: « than for x m L, i.e., pa > Pb, which result in 
fluid flow at the interface, from the right to the left. When pa 
is comparable to the nominal squeezing pressure /n/L (where 
/n is the squeezing force per unit length of the cylinder) , the 
elastic solid wall will deform and tilt as indicated in the figure. 



where we now denote L = L^. Eq. (32) agree with the 
result presented in Ref. pl| . 

For elasticahy soft materials like rubber the calculation 
of the leak-rate presented above can be simplified because 
a small change in the interfacial separation will have a 
negligible influence on the nominal stress distribution in 
the nominal contact area. Thus we can consider po{x) as 
a given fixed function obtained by squeezing the elastic 
solid against a flat surface in the absence of the fiuid. 
Using (27) and (28) we get 



Pcont(a;) = Po{x) - Pa + (Pa - Ph) 



J^dx' CF'siPcontix')) 

(33) 

This equation can be iterated to obtain the solution 
Pcont{x). If the interfacial separation u > ih-c-aia we can 
obtain the interfacial separation u from Pcont using 

U = Uo\og{l3E* /pcont), 

but in general the relation between u and Pcont must be 
calculated from (A4). 

7. Experimental 

We have performed a very simple experiment to test 
the theory presented in Sec. 6. In Fig. [TT] we show our 
set-up for measuring the leak-rate of seals. A glass (or 
PMMA) cylinder with a rubber ring (with rectangular 
cross-section) attached to one end is squeezed against a 
hard substrate with well-defined surface roughness. The 
cylinder is filled with water, and the leak-rate of the fluid 
at the rubber-countersurface is detected by the change in 
the height of the fluid in the cylinder. In this case the 
pressure difference Ap — p^ — Ph = pgH, where g is the 



glass 
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FIG. 11: Experimental set-up for measuring the leak-rate 
of seals. A glass (or PMMA) cylinder with a rubber ring 
attached to one end is squeezed against a hard substrate 
with well-defined surface roughness. The cylinder is filled 
with water, and the leak-rate of the water at the rubber- 
countersurface is detected by the change in the height of the 
water in the cylinder. 



gravitation constant, p the fluid density and H the height 
of the fluid column. With i7 w 1 m we get typically 
Ap 0.01 MPa. In the present study we use a rubber 
ring with the Young's elastic modulus E = 2.3 MPa, and 
with the inner and outer diameter 4 cm and 6 cm, re- 
spectively, and the height 0.5 cm. The rubber ring was 
made from a silicon elastomer (PDMS) prepared using 
a two-component kit (Sylgard 184) purchased from Dow 
Corning (Midland, MI). The kit consist of a base (vinyl- 
terminated polydimethylsiloxane) and a curing agent 
(methylhydrosiloxane-dimethylsiloxane copolymer) with 
a suitable catalyst. From these two components we pre- 
pared a mixture 10:1 (base/cross linker) in weight. The 
mixture was degassed to remove the trapped air induced 
by stirring from the mixing process and then poured into 
casts. The bottom of these casts was made from glass to 
obtain smooth surfaces. The samples were cured in an 
oven at 80°C for 12 h. 

We have used a sand-blasted PMMA as substrate. The 
root-mean-square roughness of the surface is 34 yum. In 
Ref. jlT| we show the height probability distribution 
P{h) and the power spectrum C{q) of the PMMA sur- 
face. 

8. Experimental results and comparison with 
theory 

In earlier studies we have performed experiments with 
the external load was so large that the condition p^ » 
Ap was satisfied, which is necessary in order to be able 
to neglect the infiuence on the contact mechanics from 
the fluid pressure at the rubber-countersurface 11 , . 




FIG. 12: Fluid leak rate as a function of the fluid pressure 
difference Ap = Pa, — Ph- The nominal squeezing pressure is 
po — 11.8 kPa. The square symbols are measured data while 
the solid lines are the theory predictions. In the calculation 
we used a = 0.54. For sandblasted PMMA with the root- 
mean-square roughness 34 fim. 

However, here we are interested in the situation where 
the fluid pressure is comparable to the nominal squeez- 
ing pressure. Thus the normal load is 18.5 N giving the 
nominal squeezing pressure po = 11.8 kPa. Using a water 
column with height H = 1.2 m gives the fluid pressure 
Pa. — Ph — 11-8 kPa at the bottom of the fluid column. 

Let us compare the theory to experiment. In Fig. [12] 
we show the fluid leak rate as a function of the fluid 
pressure difference Ap — Pa,—Ph- The square symbols are 
measured data while the solid lines are the theory predic- 
tions. Note that the fluid leak rate rapidly increases when 
the fluid pressure Ap approaches the nominal squeezing 
pressure Po = 11.8 kPa. Both the critical junction and ef- 
fective medium theories predict nearly the same pressure 
dependence of the leak-rate as observed in the experi- 
ment. 

In Fig. [13] we show the calculated contact pres- 
sure as a function of the distance x between the high- 
pressure and low-pressure side. Calculations are shown 
for (pa — Pb)/po = 0, 0.5 and 1. In Fig. [T3]we show the 
interfacial separation as a function of the distance x be- 
tween the high-pressure and low-pressure side. In Fig. [15] 
we show for p^—Ph = Po the deformed rubber block. The 
dashed line indicate the rubber block when Pa — Pb = 0. 
In this case the average separation is determined by the 
substrate surface roughness. 

9. Summary and conclusion 

In this paper we have studied fluid squeeze-out from 
the interface between an elastic solid with a flat sur- 
face and a randomly rough surface of a rigid solid. We 
have presented a very general formalism for calculating 
the (average) interfacial separation as a function of time. 
In the theory enters the effective flow conductivity cTos. 



FIG. 13: The contact pressure as a function of the distance x 
between the high-pressure and low-pressure side. The nom- 
inal squeezing pressure is 11.8 kPa. Calculations are shown 
for (pa -Pb)/po = 0, 0.5 and 1. For sandblasted PMMA with 
the root-mean-square roughness 34 ^m. 




(cm) 



FIG. 14: The interfacial separation as a function of the dis- 
tance X between the high-pressure and low-pressure side. The 
nominal squeezing pressure is po = 11-8 kPa. Calculations are 
shown for (pa— Pb)/po ~ 0, 0.5 and 1. For sandblasted PMMA 
with the root-mean-square roughness 34 /im. 



This quantity is a function of the (local) contact pressure 
Pcont- In this paper we have calculated acs using the so 
called average-separation and critical-junction theories. 
An even more accurate method, based on the Brugge- 
man effective medium theory, was developed in Ref . , 
but this theory gives results very similar to the critical- 
junction theory. The critical-junction theory and the ef- 
fective medium theory both consider the flow of fluid in 
interfacial channels and the possibility (at high enough 
squeezing pressures) of trapped fluid at the interface as 
a result of percolation of the contact area, resulting in 
confined regions (islands) of non-contact area filled with 
fluid. We have shown how this affect the time depen- 
dence of the interfacial separation. We have shown how 
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FIG. 15: Shape of the rubber block for Pa ~ Ph ~ Po- The 
nominal squeezing pressure is po ~ 11-8 kPa. The dashed line 
is the shape of the block when Pa. ^ Ph ~ For sandblasted 
PMMA with the root-mean-square roughness 34 fim. 



magnification C, 



elastic sotid 




FIG. 16: An asperity contact region observed at the magnifi- 
cation (. It appears that complete contact occur in the asper- 
ity contact region, but when the magnification is increasing 
to the highest (atomic scale) magnification i^i, it is observed 
that the solids are actually separated by the average distance 



the present theory can be used to calculate the leak-rate 
of static seals when including the reduction in the con- 
tact pressure resulting from the fluid pressure acting on 
the solids in the interfacial region. We have presented 
new experimental data which agree well with the theory 
prediction. 
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Appendix A 

Consider the elastic contact between two solids with 
randomly rough surfaces. The (apparent) relative con- 
tact area A{Q/Aq at the magnification ( can be ob- 
tained using the contact mechanics formalism developed 
elsewhere[10. Ell-Esll . where the system is studied at dif- 
ferent magnifications C. We have 13, IS] 



A{C) 1 



Ao (7rG)i/2 7q 
where 



da e--^^/^^ = erf 



2G1 



_ / p \ 2 j-Cqo 



4 \ 



where the surface roughness power spectrum 
1 



Ciq) 



(2^) 



d^x-(/i(x)/i(0))e- 



iq-x 



where (...) stands for ensemble average. Here E and v 
are the Young's elastic modulus and the Poisson ratio of 



the rubber. The height profile h{x) of the rough surface 
can be measured routinely today on all relevant length 
scales using optical and stylus experiments. 

We define Ui(C) to be the (average) height separating 
the surfaces which appear to come into contact when the 
magnification decreases from ( to — A^, where AC is a 
small (infinitesimal) change in the magnification. mi(C) 
is a monotonically decreasing function of C, and can be 
calculated from the average interfacial separation u{Q 
and ^(C) using (see Ref. _23]) 



u,{o = u{(:) + u\c)A{C)/A\c)- 



(AS) 



The quantity u{C) is the average separation between the 
surfaces in the apparent contact regions observed at the 
magnification C, see Fig. [161 It can be calculated from[23l| 

"(C) = A r dq q^C{q)w{qX) 



JpiC) P' 



where p(C) = poAq/A{C,) (where po — Pcont denote the 
nominal contact pressure) and 



w{q,C)^(nJ' dq' q''C{q')^ 



-1/2 



We will now show that as po ^ 0, for the values of 
the magnification C which are most important for the 
fiuid fiow between the solids, ui(C) — u{C,). This result 
is physically plausible because at low contact pressures 
the separation between the walls is large and the surface 
roughness should have a very small infiuence on the fluid 



flow, which therefore can be accurately studied using the 
(average) interfacial separation u for ^ « 1. 

Most of the fluid flow occur in the flow channels which 
appear close to the percolation limit where A(C)/j4o ~ 
0.5, or, using (Al), for pol G'^/'^{C,) « 1. Thus as po ^ 
we must have G[C) — >■ which, using (A2), implies C 1- 
In fact, using G{C,) ^ and (A2) one can easily show 
that for C close to unity C ~ 1 ~ Po- From (A4) it is easy 
to show that as po "> 0, the average separation will 
diverge as ^ — logpo while u'(C) diverge as —p'{C,)/p{0 — 
A'{Q)/A{C)- Thus the product u' {C,)A' {(,) / A{C,) remains 
constant as po ^ 0. It follows from (A3) that as po ^ 0, 

ui{0 ^ uiO. 

Note that from (Al) and (A2) one can calculate 

A(C)^4(C-1) ^° 

so the slope of the curve A{C,) in the relevant C-region 
becomes very high as po — ^ 0. 
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